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The identification of multiple harmonic sources (MHS) is vital to identify the 
root causes and the mitigation technique for a harmonic disturbance. This 
paper introduces an identification technique of MHS in a power distribution 
system by using a time-frequency distribution (TFD) analysis known as a 
spectrogram. The speetrogram has advantages in term of its accuracy, a less 
complex algorithm, and use of low memory size compared to previous 
methods sueh as probabilistic and harmonic power flow direetion. The 
identification of MHS is based on the significant relationship of speetral 
impedances, which are the fundamental impedance (Zi) and harmonic 
impedance (Zh) that estimate the time-frequency representation (TFR). To 
verify the performance of the proposed method, an IEEE test feeder with 
several different harmonic producing loads is simulated. It is shown that the 
suggested method is excellent with 100% correet identification of MHS. The 
method is aecurate, fast and cost-efficient in the identifieation of MHS in 
power distribution arrangement. 
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1. INTRODUCTION 

Nowadays in an advanced power distribution system, multiple harmonic sources (MHS) comprise 
nonlinear loads lead to the possible threat of a standard performance by the power system as per explained in 
[l]-[3]. Thus, it is vital to find the root cause, eliminate the harmonic pollution, and identify MHS through 
the best approach towards raising the power quality of the system [4], [5]. In addition, referring to the point 
of common coupling (PCC), MHS can be distinguished at the utility and client side which is known as 
upstream and downstream [6]-[8]. 

Various methodologies of MHS identification have been highlighted in the previous literature. 
However, there are no standards on the method to identify the dominant harmonic disturbances [30]. The 
earliest method to identify MHS is by using a power direction technique. However, this technique could not 
produce precise outcomes when the harmonics and phase angle were beyond 90° [9-10]. Furthermore, a state 
estimation method is introduced with the execution of least square and Kirchhoff current law [11]. 
Nevertheless, a study [12] depicted that state estimation technique demands many measurement devices and 
a high setup cost for a large power system. Moreover, in the literature, the identification of MHS using the 
critical impedance method has also been discussed [13-16]. The crucial drawback of this technique is a pre¬ 
requirement of source internal impedance values on the utility and client side. In practice, it is difficult to get 
those parameters without switching tests, which is costly [19-20]. To overcome the restriction of the critical 
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impedance method, a harmonic current vector method is presented [21-23]. The measurement of client 
impedances can be performed without any switching testing. However, it needs long-term measurement and 
it is difficult to analyse harmonic impedances at the client side [10]. In earlier studies [24-25], the harmonic 
source estimation has been suggested by employing the harmonic state estimation (HSE) and Bayesian 
estimation. Unfortunately, this method requires a high multiplicity in the algorithmic and requires a high cost 
to set up the distributed measurement system station [26]. By using an alternative methodology that 
concentrates on two-point measurement strategy, impedances data can be obtained by comparing the data 
between the supply and incoming load utilizing frequency and phase approach [16-18]. This approach gives 
accurate impedance data, but then again, it is still high cost and difficult to practice [20]. In reality, accuracy, 
quick estimation and low cost are important elements to consider. However, limited studies have focused on 
these factors. Based on this discussion, the restriction of previous techniques can be overcome by employing 
time-frequency distribution (TED). 

This paper proposes a high accuracy, fast estimation and cost-efficient method to identify MHS 
utilizing a single-point measurement [8] at the PCC, utilizing a TED which is spectrogram. The spectrogram 
is a suitable method that presents a signal in time and frequency representation which is known as time- 
frequency representation (TER) [9], [16], [20], [30]. Erom these signals, spectral data can be monitored 
through changes of time. Eurthermore, spectral information of the system can be extracted from TER to 
identify the system characteristics and a common harmonic producing load which is an inverter will be used 
as a load in the system [3-4]. Erom the TER, spectral impedances, which are fundamental impedance Zy, and 
harmonic impedance Z/^, will be calculated accurately. Einally, utilizing the significant relationship between 
Zy and Z/j, the identification of MHS can be distinguished in a fast and accurate manner. 


2. METHODOLOGY 

A single point-measurement of voltage and current is implemented in this research, while 
spectrogram TED is utilized for analyzing the signal parameters. 

2.1. Spectrogram 

The spectrogram is one of the techniques in TED analysis and represents a dispersion of signal 
energy in time-frequency representation. It is a squared magnitude of STET [31] thus an improvement of EET 
to cater non-stationary signals [29]. By windowing the signal at first and then taking the Eourier transforms, 
the time-localization can be obtained correctly. Eurthermore, this will lead to the rise of short time Eourier 
transform, (STET) or windowed Eourier transform. Eurthermore, the squared value of the STET is usually 
used in the signal analysis, which is called a spectrogram and can be defined as in Equation 1 and Equation 2 
[ 20 ]. 


Sy(t,f)=\J^^v(T)w(T - t)e dt\^ (1) 

tWw(T - dtp (2) 

where w(t) is observation windows. The Hanning window is selected with respect to lEC standards that 
characterize this window as suitable for analyzing harmonics [3] and is also chosen due to its lower peak side 
lope which has a narrow effect on other frequencies around fundamental value and other frequency 
components. The spectrogram is used in numerous applications and utilized as an initial investigative 
instrument as it has the property of non-producing cross-term parameters of TED [5], [6]. 

2.2. Implementation of Proposed Technique 

The implementation for MHS identification can be realized as demonstrated in Eigure 1 and Eigure 
2. Referring to the PCC, there are four cases considered in this study as follows: 

• No MHS within the distribution system. 

• Downstream: When MHS located at customer side. 

• Upstream: When MHS located at utility side. 

• Both streams: When MHS located between utility and customer side. 

Eirst of all, the measurement data are analyzed using spectrogram and presented TER. Eurthermore, 
the parameters of the power system will be estimated from the power spectrum, and finally, the identification 
of harmonic sources is based on the significant relationship of spectral impedance, which is Z1 and Zh at the 
harmonic frequency. Meanwhile, to validate the proposed method, two different converter types, which are 
an AC-DC-AC PWM and a three-phase two-level single bridge converter with sinusoidal PWM modulation 
were used as a harmonic producing load (HPL) alternately in an IEEE 4-bus test feeder. 
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Figure 1. Case 1 Figure 2. IEEE 4-bus test feeder for Case 2, 3 and 4 


Whereby N is non-HPL and H is HPL 

The implementation of the methodology in executing the proposed approach can be seen in Figure 
3. From the proposed technique, it is demonstrated that the methodology can be effortlessly actualized in the 
system. This is due to the fact that there is no requirement for any complex calculation of signal parameters. 



Figure 3. The implementation procedure in the identification of multiple harmonic sources at the PCC 


2.3. Signal Parameters 

Parameters of power quality signals are estimated from the TER in order to provide the information 
of the signal in time. Meanwhile, the minimum width of the measurement window is ten cycles [27]. 


2.3.1. Instantaneous RMS Voltage 

The RMS voltage, Vrms^ can also be derived from the TER in time and referred to as instantaneous 
RMS voltage, VrmsW and written as [28], 





Pxit,f)df 


( 3 ) 


where Px(t, f) is the TER of the signal and fs is the sampling frequency of the system. 

2.3.2. Instantaneous RMS Fundamental Voltage 

Instantaneous RMS fundamental voltage, Virms(t)^ is defined as the RMS voltage at power system 
frequency [28]. From the TFR, the Vijmsit)’ can be calculated as: 


^Irms (^) 



fm = /o + 25 Hz 


( 4 ) 


^=/o-25Hz 
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Where Px(t,f) is the time-frequency representation of the signal and fO is the fundamental frequency that 
corresponds to the power system frequency. 25 Hz is chosen forf^i and//^ since it can cover the fundamental 
frequency component to calculate the magnitude of the frequency component. 

2.3.3. Instantaneous Total Harmonic Distortion 

The instantaneous total harmonic distortion, THD(t), is used to measure how much harmonic 
content is in a waveform and expressed as [28], 


h.rmsiPP 

^ irms(t) 


THD{t) = 


where Vh,rms fO is RMS harmonic voltage from 2nd to 50th harmonic. 

2.3.4. Instantaneous Total Nonharmonic Distortion 

Besides harmonic, a signal also contains interharmonic components that are not multiple integers 
and can be quantified by using the instantaneous total nonharmonic distortion, TnHD(t), an index which is 
defined as [28], 


TnHD{t) = 




where Vh,rms fO is RMS harmonic voltage from 2nd to 50th harmonic. 


( 6 ) 


2.3.5. Spectral Impedance 

The spectral impedance is defined as the ratio of the voltage phasor to the current phasor. Thus, it 
can be calculated using [12]: 


^ v-,{t)^evx 

(7) 

hit)AOn 

^ Vhit)Aevh 

(8) 

ZJt = - 

ihit)Aein 


where, /z is a harmonic number at 3, 5, 7, 9 ...n 


3. RESULTS AND DISCUSSION 

The accuracy and validity of this methodology are examined and verified in the distribution system 
with two different types of inverters. The significant relationship of spectral impedances in all cases will be 
observed due to distinguishing the location of MHS. 

3.1. Case 1: No Harmonic Sources within the System 

Eigure 4(a) demonstrates the measured voltage and current in the time domain at PCC N-N. 
Meanwhile, Eigure 4(b) shows that THD(t) and TnHD(t) are zero percent. This result demonstrates that no 
harmonic and interharmonic components exist in the signal. Eurthermore, TERs in Eigure 4(c) and 4(d) 
clearly show that only the fundamental component of voltage and current exist at 50 Hz. Erom the TER, the 
Z] is estimated with a value of 20 ohm and presented in Eigure 4(d). 

Eigure 4(e) clearly indicates the significant relationship between Zj and Zh^ that been estimated from 
TER of spectral impedance for Case 1 can be concluded as: 


Zh=0 ohm 

(9) 

ZjfO ohm 

(10) 
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Figure 4. No harmonic sources, (a) Voltage and current signal, (b) THD(t) and TnHD(t) for voltage, (c) 
TFR of voltage, (d) TFR of current, (e) Spectral impedance component 


3.2. Case 2: Harmonic Sources at Downstream 

In case 2, the effectiveness of the proposed method is analyzed when HPL is located at downstream 
and the measurement data are taken at the PCC N-H. Figure 5(a) shows the harmonic signal in the time 
domain when the AC-DC-AC PWM converter was used as HPL at downstream. Furthermore, Figure 5(b) 
illustrates that the THD(t) is in the range of 2.5% to 9.8% and it indicates that the signal contains the 
harmonic component. Meanwhile, the interharmonic index, which is TnHD (r) is in the range of 27% to 37%. 
Figure 5(c) and 5(d) show the signals of five frequency components consist of the fundamental frequency 
(50Hz), the harmonic component at 375Hz, and interharmonic components at 262.5Hz, 587.5Hz and 
687.5Hz, respectively. Figure 5 (e) presents the TFR of spectral impedance and as depicted in Figure 5(f), the 
estimated impedance value of Z; is 8 ohm and each Zh has a value of 6 ohm. It is somewhat surprising that 
each value of Zh is smaller than Z; in this case. It is possible to hypothesize that there is a possible significant 
relationship between Zy and Zh in identifying the location of the harmonic source. Therefore, to conclude the 
findings in this case, further analysis was conducted by replacing the HPL with a three-phase two-level single 
bridge converter. 
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Eigure 5. AC-DC-AC PWM converter as a harmonic source at upstream, (a) Voltage and current signal, 
(b) THD(t) and TnHD(t) for voltage (c) TER of voltage, (d) TER of current, (e) TER of spectral 
impedance, (f) Spectral impedance component 


The voltage and current signal as shown in Eigure 6 (a) is obtained when a three-phase two-level 
single bridge converter is used as HPL at downstream. Eurthermore, 6(b) illustrates that the THD(t) is in the 
range of 2 % to 8% and it shows that the signal comprises a harmonic component. Meanwhile, the TnHD(t) is 
in the range of 27% to 38%. In addition, Eigure 6(c) and 6(d) depict the TER of voltage and current. It can be 
seen from the figures that the fundamental component has 50Hz, harmonic components at 375Hz and 
interharmonic components discovered at 262.5Hz and 587.5Hz, respectively. Eurthermore, Eigure 6(e) 
presents the TER of spectral impedance. Eigure 6(f) shows the estimated value of Z; is 11 ohm, while each Zh 
has a value of 8.5 ohm. What is surprising is that each value of Zh is smaller than Z; which is similar to the 
previous result. The outcomes of this study have proven the significant relationship of spectral impedance in 
identifying the MHS. 
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Figure. 6. Three-phase two-level single bridge converter as a harmonic source at upstream, (a) Voltage 
and current signal, (b) THD(t) and TnHD(t) for voltage, (c) TFR of voltage, (d) TFR of current, (e) TFR 
of spectral impedance, (f) Spectral impedance component 


The results obtained from Figure 5(f) and Figure 6(f) of Case 2 clearly indicates that Z; value is 
always greater than Z^. Therefore, the significant relationship of spectral impedance for Case 2 can be 
concluded as: 

Zi>Zn (11) 

3.3. Case 3: Harmonic Sources at Upstream and Downstream 

The goal of this case study is to localize the harmonic sources when both streams contain the HPL 
and the measurement is taken at PCC H-H. It is hypothesized that another significant relationship of spectral 
impedance will be identified from this case. Current and voltage signals of this case are shown in Figure 7(a) 
when an AC-DC-AC PWM converter acted as HPL at upstream and downstream. As can be seen from Figure 
7(b), the THD(t) is in the range of 3% to 10% and it indicates that the signal contains the harmonic 
component. Meanwhile, the interharmonic index which is TnHD (t) is in the range of 27% to 37.5%. Figure 
7(c) and 7(d) show the signals of five frequency components consist of the fundamental frequency (50Hz), 
the harmonic component at 375Hz, and interharmonic components at 262.5Hz, 587.5Hz and 687.5Hz, 
respectively. In addition. Figure 7(e) presents the TFR of spectral impedance. As shown in Figure 7(f), the 
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estimated spectral impedance from TER comprise of Z; is 13 ohm, while each Z/^ has a value of 14.5 ohm. 
Interestingly, is it observed that each value of is greater than Z; when the MHS are located at downstream 
and upstream. To strengthen the preliminary findings, the MHS was replaced with a three-phase two-level 
single bridge converter and properly analyzed in the next discussion. 
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Eigure 7. AC-DC-AC PWM converter as a harmonic source at upstream and downstream, (a) Voltage and 
current signal, (b) THD(t) and TnHD(t) for voltage, (c) TER of voltage, (d) TER of current, (e) TER of 
spectral impedance, (f) Spectral impedance component 


Eor further investigation of Case 3, the three-phase two-level single bridge converter is employed as 
HPL at upstream and downstream. Eigure 8(a) shows the signal of voltage and current in the time domain at 
PCC H-H. Moreover, Eigure 8(b) shows that the THD(t) is in the range of 2.5% to 8% and it indicates that 
the signal contains the harmonic component. Meanwhile, the interharmonic index, shown as TnHD (t) is in 
the range of 21% to 31%. Eigure 8(c) and 8(d) show the signals comprise a fundamental frequency at 50Hz, 
the harmonic component at 375Hz, and interharmonic components at 262.5Hz and 587.5Hz, respectively. 
Besides, Eigure 8(e) depicts the TER of the spectral impedance, and it is apparent from Eigure 8(f) that the Zj 
is estimated with a value of 18 ohm and each Zh has a value of 22 ohm. Interestingly, the spectral impedance 
relationship, in this case, is similar even when using a different type of converter such as the HPL, whereby 
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each of Z/^ is greater than the Z 7 . It is interesting to observe that in this case, the significant spectral 
impedance relationship can be obtained from the proposed method and utilized in identifying the MHS. 
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Figure 8 . Three-phase two-level single bridge converter as a harmonic source at upstream and downstream, 
(a) Voltage and current signal, (b) THD(t) and TnHD(t) for voltage, (c) TFR of voltage, (d) TFR of current, 
(e) TFR of spectral impedance, (f) Spectral impedance component 


The evidence from Figure 7(f) and Figure 8 (f) clearly indicate that the Z1 is always smaller than the 
Zh when the MHS is located at downstream and upstream. Thus, the significant relationship of spectral 
impedances for Case 4 can be written as: 

Zj<Zh ( 12 ) 

3.4. Case 4: Harmonic Sources at Upstream 

With the aim of identifying the HPL located upstream. Case 4 is executed with an AC-DC-AC 
PWM converter as HPL downstream, and a measurement is taken at PCC H-N. Figure 9(a) shows the voltage 
and current signal in time domain. As illustrated in Figure 9(b), the THD(t) is in the range of 2.5% to 10%, 
indicating that the signal contains the harmonic component. Meanwhile, the range of TnHD (t) is between 
28% to 38%. In addition. Figure 9(c) and 9(d) show the components of the signal that consist of the 
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fundamental frequency (50Hz), the harmonic component at 375Hz, and interharmonic components at 
262.5Hz, 587.5Hz and 687.5Hz, respectively. Eigure 9(e) also presents the TER of spectral impedance and 
Eigure 9(f) shows that the estimated value of Z1 and Zh are equal with a value of 40 ohm. It is, therefore, 
likely that there is a relationship of spectral impedance in this case. To verify this hypothesis, the HPL was 
replaced with a three-phase two-level single bridge converter. 
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Eigure 9. AC-DC-AC PWM converter as a harmonic source at downstream, (a) Voltage and current signal, 
(b) THD(t) and TnHD(t) for voltage, (c) TER of voltage, (d) TER of current, (e) TER of spectral 
impedance, (f) Spectral impedance component 


A three-phase two-level single bridge converter replaces the previous MHS and the time domain 
signals are as shown in Eigure 10(a). In the meantime, Eigure 10(b) shows that the THD(t) is in the range of 
2.5% to 10% and the TnHD (t) is between 27% to 37%. Eigure 10(c) and 10(d) show the signals consist of a 
fundamental frequency of 50Hz, the harmonic component at 375Hz, and interharmonic components at 
262.5Hz and 587.5Hz, respectively. Eigure 10(e) also depicts the TER of the spectral impedance. As can be 
observed in Eigure 10(f), the Zj and Zh were estimated at 40 ohm. The results indicate that the value of the Zj 
and Zh are constantly equal in this case. 
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Figure 10. Three-phase two-level single bridge converter as a harmonic source at downstream, (a) Voltage 
and current signal, (b) THD(t) and TnHD(t) for voltage, (c) TFR of voltage, (d) TFR of current, (e) TFR 
of spectral impedance, (f) Spectral impedance component 


It is apparent from Figure 9(f) and Figure 10(f), that the value of the Zj and Zh are constantly similar 
for Case 4. Therefore, the significant relationship of spectral impedances for Case 4 can be concluded as: 

Z;=Z, (13) 


Table 1. The S ignificant Relationship of Spectral Impedance for MH S Localisation 


Case 

Upstream 

Downstream 

Spectral Impedance (Q) 

1 

N 

N 

Zh=0 

Zj^ 

2 

N 

H 

Zj>Zh 

3 

H 

H 

Zi<Zh 

4 

H 

N 

Zi=Zh 


Where by N is non-HPL, H is HPL, the Zj is fundamental impedance and Zh is harmonic impedance. 
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4. CONCLUSION 

The current findings in this study indicate the significant relationship and noteworthy contributions 
of spectral impedance in MHS identification. The previous methods such as power direction method and 
probabilistic method that require a high cost to set up a multi-point measurement station with high 
multiplicity use of the algorithm. However, the proposed method only requires a single point-measurement 
that utilizes a fast and simple algorithm that is short time Eourier transform-based. 

The validation of the proposed method is done by using two different types of MHS. Based on Table 
1, MHS in the power system can be identified and distinguished clearly in a significant pattern as follows: 

• If the Z/j is 0 ohm, there is no MHS within the system. 

• If the Zy is greater than Z/^, there is MHS at downstream. 

• If the Zy is smaller than Z/^, there are MHS at upstream and downstream. 

• If the Zy is equal to Z/^, there is MHS at upstream. 

Thus, the improved methodology proposed can be used accurately and quickly to obtain the 
significant relationship of spectral impedances in order to distinguish the location of MHS. In addition, the 
proposed method is excellent and provides 100% correct identification of MHS in the power system. 
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